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No Evidence for Dark Energy Dynamics from 
a Global Analysis of Cosmological Data 

We use a variant of principal component analysis to investigate the possible temporal evolution 
\^ . of the dark energy equation of state, w(z). We constrain w(z) in multiple redshift bins, 

pH ' utilizing the most recent data from Type la supernovae, the cosmic microwave background, 

baryon acoustic oscillations, the integrated Sachs- Wolfe effect, galaxy clustering, and weak 
lensing data. Unlike other recent analyses, we find no significant evidence for evolving dark 
energy; the data remains completely consistent with a cosmological constant. We also study 
the extent to which the time-evolution of the equation of state would be constrained by a 
combination of current and future-generation surveys, such as Planck and the Joint Dark 
Energy Mission. 



1 Introduction 



■^j- ■ One of the denning challenges for modern cosmology is understanding the physical mechanism 

responsible for the accelerating expansion of the Universe > 2 . The origin of the cosmic acceler- 
ation can be due to a new source of stress-energy, called "dark energy", a modified theory of 
gravity, or some mixture of both 3 . 

In the absence of a well-defined and theoretically motivated model for dark energy, it is gen- 
erally assumed that the dark energy equation of state (the ratio of pressure to energy density) 
evolves with redshift with an arbitrary functional form. Common parameterizations include a 
linear variation, w{z) = wq + w z z A , or an evolution that asymptotes to a constant w at high 
redshift, w(z) = wq + w a z/(l + z) 5 ' 6 . However, given our complete ignorance of the underlying 
physical processes, it is advisable to approach our analysis of dark energy with a minimum of 
assumptions. Fixing an ad hoc two parameter form could lead to bias in our inference of the 
dark energy properties. 

In this paper we measure the evolution history of the dark energy using a flexible and almost 
completely model independent approach, based on a variant of the principal component analysis 
(PCA) introduced in Huterer (2003) 7 ; in order to be conservative, we begin by using data we 
determine the equation only from geometric probes of dark energy, namely the cosmic microwave 
background radiation (CMB), Type la supernovae (SNe) and baryon acoustic oscillation data 
(BAO). We perform a full likelihood analysis using the Markov Chain Monte Carlo approach 8 . 
We then consider constraints on w(z) from a larger combination of datasets, including probes of 
the growth of cosmological perturbations, such as large scale structure (LSS) data. An important 
consideration for such an analysis is to properly take into account dark energy perturbations, 
and we make use of the prescription introduced in 9 . This method implements a Parameterized 
Post-Friedmann (PPF) prescription for the dark energy perturbations following 10 . 



2 Analysis and results 

The method we use to constrain the dark energy evolution is based on a modified version of 
the publicly available Markov Chain Monte Carlo package CosmoMC 8 , with a convergence 
diagnostics based on the Gelman-Rubin criterion 12 . We consider a flat cosmological model 
described by the following set of parameters: 

{wi,uj b ,uj c ,@s,T,n s ,log[l0 10 A s ]} , (1) 

where u^ (= f^/i 2 ) and ui c (= £l c h 2 ) are the physical baryon and cold dark matter densities 
relative to the critical density, Q s is the ratio of the sound horizon to the angular diameter 
distance at decoupling, r is the optical depth to re-ionization, and A s and n s are the amplitude 
of the primordial spectrum and the spectral index, respectively. We bin the dark energy equation 
of state in five redshift bins, Wi(z) (i = 1,2, ..5), representing the value at five redshifts, Z{ G 
[0.0, 0.25, 0.50, 0.75, 1.0] and, for z > 1, we fix the equation of state parameter at its z = 1 value, 
since we find that current data place only weak constraints on w(z) for z > 1. To summarize, 
our parameterization is given by: 



w(z) = { Wi, z <z max ,z € {zi}; (2) 



w{z = 1), 


z > 1; 


Wi, 


Z _; Zmax j Z fc \ Zi j , 


spline, 


Z _; Zmax j Z (jt \ Zi j . 



Finally, we follow n to determine uncorrelated estimates of the dark energy parameters. 
For the CMB, we use data and likelihood code from the WMAP team's 5-year release 13 (both 
temperature TT and polarization TE; we will refer to this analysis as WMAP5). We also checked 
that results don't change if we use the latest data release from WMAP 14 . Supernova data come 
from the Union data set (UNION) produced by the Supernova Cosmology Project 15 ; however, 
to check the consistency of our results, we also used the recently released Constitution dataset 
(Constitution) 18 which, with 397 Type la supernovae, is the largest sample to date. We also 
used the latest SDSS release (DR7) BAO distance scale 16 ' 17 . Weak lensing (WL) data are taken 
from CFHTLS 22 and we use the weak lensing module provided in 19 > 20 ; with some modifications 
to assess the likelihood in terms of the variance of the aperture mass (Eq. 5 of 22 ) with the full 
covariance matrix 21 . The cross-correlation between CMB and galaxy survey data is employed 
using the public code at 23 . We modify it to take into account the temporal evolution of the dark 
energy equation of state, since the code only considers wCDM cosmologies. We refer to 24 ' 25 for 
a description of both the methodology and the datasets used. Finally, we use the recent value of 
the Hubble constant from the SHOES (Supernovae and Hq for the Equation of State) program, 
Hq = 74.2 ± 3.6 km s _1 Mpc~ 1 (l<r) 26 . We also incorporate baryon density information from 
Big Bang Nucleosynthesis Qbh 2 = 0.022 ± 0.002 (la) , as well as a top-hat prior on the age of 
the Universe, 10 Gyr < t < 20 Gyr. 

As we can see from Fig. 1, all values are compatible with a cosmological constant (w = — 1) at 
the 2<r level; in particular, there is no discrepancy between the Union and Constitution datasets. 
Moreover, as we can see from Fig. 2, the addition of cosmological probes of cosmic clustering 
noticeably reduces the uncertainty in the determination of the dark energy parameters, especially 
at high redshifts. 

To reinforce our conclusions, we also created several mock datasets for upcoming and future 
SN, BAO, and CMB experiments. The quality of future datasets allows us to constrain the 
dark energy evolution beyond redshift z = 1. We thus consider an additional bin at z = 1.7, 
with a similar constraint: w(z > 1.7) = w(z = 1.7). We consider a mock catalog of 2,298 SNe, 
with 300 SNe uniformly distributed out to z = 0.1, as expected from ground-based low redshift 
samples, and an additional 1998 SNe binned in 32 redshift bins in the range 0.1 < z < 1.7, as 
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Figure 1: Left: uncorrelated constraints on the dark equation of state parameters using WMAP+UNION+BAO. 
Right: comparison between WMAP+UNION+BAO and WMAP+Constitution+BAO; the points for the Consti- 
tution dataset have been slightly shifted to facilitate comparison between the two cases: we find no significant 

difference between UNION and Constitution. 
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Figure 2: Results using data from a "global" dataset which includes WMAP+UNION+BAO+WL+ISW+LSS; 

error bars are at 2a. 



expected from JDEM or similar future surveys 28 . The error in the distance modulus for each 
SN is given by the intrinsic error, cr- m ^ = 0.1 mag. In addition, we use a mock catalog of 13 BAO 
estimates, including 2 BAO estimates at z = 0.2 and z = 0.35, with 6% and 4.7% uncertainties (in 
Dy), respectively, 4 BAO constraints at z = [0.6,0.8, 1.0, 1.2] with corresponding fiducial survey 
precisions (in Dy) of [1.9,1.5,1.0,0.9]% (V5N5 from 29 ), and 7 BAO estimates with precision 
[0.36, 0.33, 0.34, 0.33, 0.31, 0.33, 0.32]% from z = 1.05 to z = 1.65 in steps of 0.1 30 . 
We simulate Planck data using a fiducial ACDM model, with the best fit parameters from 
WMAP5, and noise properties consistent with a combination of the Planck 100-143-217 GHz 
channels of the HFI , and fitting for temperature and polarization using the full-sky likelihood 
function given in 32 . In addition, we use the same priors on the Hubble parameter and on the 
baryon density as considered above. As can be seen from Table 1 and Figure 3, future data will 
reduce the uncertainties in Wi by a factor of at least 2, with the relative uncertainty below 10% 
in all but the last bin (at z = 1.7). 
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Figure 3: Uncorrelated constraints on the dark energy equation of state parameters, for mock datasets from 

Planck and JDEM; error bars are at 2<r. 

3 Conclusions 

One of the main tasks for present and future dark energy surveys is to determine whether or 
not the dark energy density is evolving with time. 

We have performed a global analysis of the latest cosmological datasets and have constrained the 
dark energy equation of state using a very flexible and almost model independent parameteriza- 
tion. We determine the equation of state w(z) in five independent redshift bins, incorporating 
the effects of dark energy perturbations. We find no evidence for a temporal evolution of dark 
energy — the data is completely consistent with a cosmological constant. This agrees with most 
previous results, but significantly improves the overall constraints. We show that future exper- 
iments, such as Planck or JDEM, will be able to reduce the uncertainty on w(z) to less than 
10% in multiple redshift bins, thereby mapping any temporal evolution of dark energy with high 
precision. With this data it will be possible to measure the temporal derivative of the equation 
of state parameters, dw/dz, useful in discriminating between two broad classes of "thawing" 
and "freezing" models 33 . 
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